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The  construcal  tree-shaped  network  is  introduced  into  the  design  of  a  methanol  steam  microreactor  in 
the  context  of  optimization  of  the  flow  configuration.  A  three-dimensional  model  for  methanol  steam 
reaction  in  this  designed  microreactor  is  developed  and  numerically  analyzed.  The  methanol  conver¬ 
sion,  CO  concentration  in  the  product  and  the  total  pressure  drop  of  the  gases  in  the  microreactor  with 
constructal  tree-shaped  network  are  evaluated  and  compared  with  those  in  the  serpentine  reactor.  It  is 
found  that  the  reaction  of  methanol  steam  reforming  is  enhanced  in  the  constructal  tree-shaped  microre¬ 
actor,  since  the  tree-shaped  reactor  configuration,  which  acts  an  optimizer  for  the  reactant  distribution, 
provides  a  reaction  space  with  larger  surface-to-volume  ratio  and  the  reduction  of  reactant  velocities 
in  the  branches.  Compared  with  the  serpentine  microreactor,  the  constructal  reactor  possesses  a  higher 
methanol  conversion  rate  accompanied  with  a  higher  CO  concentration.  The  conversion  rate  of  the  con¬ 
structal  microreactor  is  more  than  10%  over  that  of  serpentine  reactor.  More  particularly,  the  reduction 
of  flow  distance  makes  the  constructal  microreactor  still  possess  almost  the  same  pressure  drop  as  the 
corresponding  serpentine  reactor,  despite  that  the  bifurcations  induce  extra  local  pressure  loss,  and  the 
reduction  of  channel  size  in  branches  also  causes  pressure  losses. 
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1.  Introduction 

There  is  an  increasing  interest  in  using  methanol  as  the  fuel  for 
fuel  cells,  whether  used  directly  [1,2]  or  preprocessed  by  steam 
reforming  [3],  partial  oxidation  [4],  or  autothermal  reforming  [5]. 
For  the  applications  of  portable  electrical  power,  of  special  interest 
is  methanol  steam  reforming  to  produce  hydrogen  for  the  proton 
exchange  membrane  fuel  cells  (PEMFC),  due  to  its  high  volumet¬ 
ric  energy  densities,  relatively  low  operating  temperatures,  high 
hydrogen  yield  and  low  CO  production  [3,6,7].  One  of  the  most 
important  challenges  for  the  use  of  PEMFC  with  methanol  steam 
reforming  in  portable  applications  is  the  availability  of  a  com¬ 
pact  and  light  unit  of  fuel  reformer  with  high  hydrogen  production 
efficiency.  Fortunately,  microreaction  technology  provided  a  good 
solution  to  meet  this  challenge  since  it  not  only  can  reduce  the  vol¬ 
ume  of  the  chemical  devices  but  also  does  significantly  enhance 
the  reaction  process  and  mass  transfer  [8].  In  particular,  a  lot  of 
methanol  steam  microreformers  have  been  designed  and  investi¬ 
gated  to  supply  hydrogen  in  the  past  decades  [3,9]. 

Kawamura  et  al.  [10]  designed  a  microreformer  with  serpen¬ 
tine  microchannel  for  the  supply  of  hydrogen  for  a  small  PEMFC. 
This  micro-reformer  demonstrated  to  be  capable  of  maintaining  a 
hydrogen  production  rate  suitable  for  powering  1  W-class  devices 
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such  as  cellular  phones.  Park  et  al.  [11]  fabricated  a  methanol 
steam  microreformer  with  serpentine  microchannel  on  a  silicon 
wafer  by  MEMS  facilities.  The  reforming  rate,  heat  transfer  and 
fluid  flow  through  the  reformer  are  experimentally  and  numerically 
investigated.  Chen  et  al.  [12]  designed  a  plate-type  micro  reformer 
constructed  by  parallel  microchannels  with  a  diagonal  inlet/outlet 
and  numerically  analyzed  the  effects  of  liquid  feed  rate,  reaction 
temperature,  and  steam  to  methanol  ratio  on  the  concentration 
profiles  of  methanol,  hydrogen,  and  carbon  monoxide  in  the  reac¬ 
tor.  Lee  et  al.  [13]  compared  the  reacting  flow  in  a  packed-bed 
reactor  and  a  wall-coated  reactor  made  of  borosilicate  glass  tubes 
with  an  inner/outer  diameter  of  1.5  mm/1.8  mm.  They  found  that 
the  wall-coated  reformer  has  a  smaller  power  requirement  for 
delivering  fuel  than  that  of  packed  catalytic  bed  reformer.  And 
the  thermal  resistance  of  coated-catalytic  layer  is  also  smaller. 
FIsueh  et  al.  [14]  numerically  investigated  a  plate  steam  methanol 
microreformer  with  parallel  microchnannels.  The  effects  of  geo¬ 
metric  sizes  (aspect  ratios  and  micro  channel  size)  and  thermo-fluid 
parameters  (Reynolds  number  and  wall  temperature)  on  methanol 
conversion  and  local  transport  phenomena  were  explored. 

Kundu  et  al.  [15]  compared  the  performance  between  a  par¬ 
allel  channeled  and  a  serpentine  channeled  microreformer.  It 
was  indicated  that  the  performance  of  the  serpentine  channeled 
micro-reformer  is  always  higher  than  that  of  parallel  channeled 
micro-reformer.  Pattekar  and  Kothare  [16]  developed  a  novel 
radial  flow  packed-bed  microreactor  for  hydrogen  generation.  The 
microreactor  can  provide  an  order  of  magnitude  reduction  in  the 
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Nomenclature 

A,  B 

Pre-exponential  term  in  the  Arrhenius  expression 

Ci 

Mole  concentration  of  species  i,  mol  m-3 

cp 

Specific  heat,  J  kg-1  K 

Cd.Cr 

Modification  factors  of  decomposition  and  reform¬ 
ing  reaction 

dk 

Branch  diameter  of  the  kth  level,  m 

D 

Length  dimension  of  channel 

£d,  £r 

Activation  energy  in  Arrhenius  expression  for 
decomposition  and  reforming  reaction,  J  mol-1 

h 

Depth  of  channels,  m 

hj 

Enthalpy  of  species  j 

AH 

Enthalpy  of  reaction,  J  mol-1 

1 

Diffusion  flux  of  species  i,  kgnrr1  K-1 

k, 

Thermal  conductivity  of  species  i,  W  m-1  K-1 

km 

Thermal  conductivity  of  the  mixture,  W  m- 1 1<-1 

kD,  kR,  ks 

Rate  constant  for  decomposition  reaction,  reform¬ 
ing  reaction  and  shift  reaction 

Kn 

Knudsen  number 

l 

Characteristic  geometric  dimension  of  the  channel, 

4 

m 

Branch  length  of  the  kth  level,  m 

f-X*  Fy,  LZ 

Length  of  the  microreactors  in  the  x,  y,  z  direction 

(m) 

M 

Molecular  weight,  kg  mol-1 

P 

Pressure,  Pa 

n 

Molar  rate  of  creation/destruction  of  species  i, 
mols-1  m2 

Rate  of  decomposition  reaction,  mol  s-1  m-2 

Jr 

Rate  of  reforming  reaction,  mol  s-1  m-2 

Re 

Reynolds  number 

SMR 

Steam  to  methanol  molar  ratio 

T 

Temperature,  K 

u 

Velocity,  ms-1 

U[n 

Inlet  velocity  of  reactants,  m  s-1 

w 

Width  of  channels,  m 

Wcat 

Catalyst  density,  kg  m-2 

X 

Conversion  of  methanol 

yi 

Mole  fraction  of  species  i 

Yi 

Mass  fraction  of  species  i 

Greek  symbols 

A 

Diameter  dimension 

X 

Molecular  mean  free  path,  m 

A i 

Viscosity,  kgnrr1  s-1 

Subscripts 

cat 

Catalyst 

D 

Decomposition 

i 

Species  index  ( 1  -CH3 OH,  2-H20, 3-H2 ;  4-CO,  5-C02 ) 

in 

Inlet 

k 

Branching  level 

m 

Mixture 

out 

Outlet 

R 

Reforming 

S 

Shift 

w 

Wall 

required  driving  pressure  while  maintaining  the  compact  design 
of  the  microreactor  and  providing  greater  than  98%  conversion  of 
methanol  to  hydrogen  at  flow  rates  sufficient  to  supply  fuel  to  a 
20  W  fuel  cell. 


Considering  that  the  flow  configuration  of  a  reactor  is  one  of  the 
most  important  factors  to  its  performance,  an  appropriate  reac¬ 
tor  flow  configuration  will  undoubtedly  improve  the  reactant  gas 
transport  and  the  efficiency  of  thermal  management,  which  results 
in  an  improvement  in  the  efficiency  of  reforming  process.  However, 
until  now,  the  most  commonly  used  flow  configurations  in  microre¬ 
actor  still  rest  on  single  channel  or  parallel  channels.  Therefore,  it  is 
desirable  to  develop  a  microreactor  configuration  with  optimized 
conversion  and  flow  performance. 

Inspiration  for  the  optimization  of  many  engineering  prob¬ 
lems  can  always  be  found  in  nature.  For  example,  advanced 
conversion  and  flow  efficiency  can  be  found  in  human  circulatory 
and  respiratory  systems.  The  flow  architectures  of  the  blood- 
circulating  arteries  and  veins  of  human  beings  are  self-organized  as 
a  branching  vessel  tree  system,  and  the  advanced  energy  and  mass 
transfer  efficiency  performs  in  these  architectures.  Constructing 
such  branch-tree  like  architectures  is  an  essential  objective  of  con- 
structal  theory  which  is  to  describe  the  geometry  and  evolution  of 
optimized  and  organized  natural  phenomena.  The  principal  of  the 
constructal  theory  was  formulated  by  Bejan  in  1 996  as  the  construc- 
tal  law  of  the  generation  of  flow  configuration:  “For  a  finite-size 
flow  system  to  persist  in  time  (to  survive)  its  configuration  must 
evolve  in  such  a  way  that  it  provides  an  easier  access  to  the  currents 
that  flow  through  it”  [17,18].  Recent  researches  have  demonstrated 
that  this  constructal  theory  [19,20]  has  been  successfully  used  to 
optimize  the  flow  configuration.  The  important  applications  of  the 
constructal  tree-shaped  architectures  have  been  found  in  electron¬ 
ics  cooling  [21-24],  fuel  cells  [25-27],  traffic  [28],  etc.  Therefore, 
should  we  construct  similar  architectures  to  optimize  the  reac¬ 
tor,  more  advanced  transport  and  reaction  performance  may  be 
acquired. 

In  the  context,  the  tree-shaped  architecture  is  introduced 
to  design  the  methanol  steam  microreactor  in  this  work.  A 
three-dimensional  model  for  methanol  steam  reaction  in  the 
microreactor  with  constructal  tree-shaped  network  is  developed 
and  numerically  analyzed.  The  methanol  conversion,  CO  concen¬ 
tration  in  the  product  and  pressure  drop  of  the  gases  in  the 
microreactor  with  constructal  tree-shaped  network  are  evaluated 
and  compared  with  those  in  the  microreactor  with  serpentine  flow 
pattern.  In  particular,  the  significant  enhancement  of  chemical 
reaction  in  tree-shaped  channels  is  identified  and  discussed. 

2.  Construction  of  a  tree-shaped  network  flow 
configuration  for  microreactors 

As  shown  in  Fig.  1,  a  rectangular  shape  tree-shaped  network 
flow  configuration  with  6  branching  levels  can  be  constructed  as 
[22,23]: 

a.  Suppose  that  every  channel  is  divided  into  two  branches  at  the 
next  level,  so  N= 2,  the  branching  angles  is  180°,  and  the  tree¬ 
shaped  network  has  6  branching  levels  as  shown  in  Fig.  1(a). 

b.  Define  that  the  ratio  of  the  length  of  the  channel  at  the  (k  + 1  )th 
branching  level  to  the  length  of  the  channel  at  the  kth  branching 
level  as 

hp±=N-vD  (i) 

H< 

It  follows  that 

Lk=L0N~k/D  (2) 

where  L0  is  the  length  at  the  0th  branching  level.  As  suggested  by 
Chen  and  Cheng  [22],  the  energy  and  mass  transfer  efficiency  at 
the  length  dimension  D  =  2  will  be  the  highest,  so  D  =  2  is  chosen 
in  this  work. 
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Fig.  1.  The  schematic  of  constructal  tree-shaped  network:  (a)  tree-shaped  network  and  (b)  branches. 


c.  The  branch  hydraulic  diameters  before  and  after  bifurcation  are 

denoted  by  dk  and  d/<+1 ,  respectively,  and  the  hydraulic  diameter 

ratio  is  defined  as 

^±1=n-VA  (3) 

dk 

where  the  diameter  dimension  A  =  3  is  chosen  according  to  the 

Murray  Law  [29].  It  follows  that 

dk  =  d0N~k/A  (4) 

where  d0  is  the  initial  hydraulic  diameter  of  channel  network. 

In  the  tree-shaped  network  configuration,  the  reduction  of 
channel  size  leads  to  a  larger  reaction  surface  in  a  specified  volume, 
which  is  beneficial  to  chemical  reaction.  However,  the  collection 
of  reactants  is  difficult  to  be  realized  in  a  single  layer  tree-shaped 
net.  To  solve  this  problem,  the  reactor  is  designed  with  two  layers 
configuration  (i.e.  inlet  and  outlet  layers  of  constructal  tree-shaped 
network),  and  the  inlet  layer  is  of  the  same  distribution  of  outlet 
layer  except  that  the  Oth  branches  in  the  inlet  and  outlet  layers  are 
in  the  opposite  direction.  The  ends  (the  highest  branching  (/<  =  6)) 
of  these  two  layers  are  communicated  by  channels  between  them 
as  shown  in  Fig.  2(a)  and  (b).  The  cross  section  of  each  branch  is 
rectangular  and  has  the  same  depth  of  h  =  1  mm. 

To  compare  the  performance  of  the  microreactor  with  con¬ 
structal  tree-shaped  network,  a  microcreactor  with  traditional 
serpentine  flow  configuration  as  shown  in  Fig.  2(c)  with  both  the 
same  channel  volume  and  overall  dimensions  (Lx,  Ly,  Lz  as  shown  in 
Fig.  2)  is  designed.  In  addition,  the  serpentine  flow  pattern  is  also 
of  the  same  cross  section  with  the  initial  channel  in  tree  network, 
i.e.  the  width  w  =  2  mm,  the  depth  h  =  1  mm  and  the  total  length  of 
the  serpentine  channel  is  110  mm.  The  geometry  parameters  for 
the  two  reactors  are  shown  in  Table  1. 

Considering  the  MEMS  technology  is  widely  and  easily  imple¬ 
mented  in  the  fabrication  of  microreactors,  the  silicon  wafer  is  used 
as  the  substrate  of  the  microreactors  in  the  simulation.  A  thin  layer 
of  commercial  catalyst  F3-01(CuO/ZnO/Al2O3)  with  a  thickness  of 
30  [xm  is  assumed  to  be  uniformly  deposited  on  channel  walls. 

3.  Methanol  steam  reforming  in  microreactor 

3.1.  Mathematical  model 

The  mixture  of  steam  and  methanol  flows  in  the  microre¬ 
actor  with  a  specified  inlet  flow  rate,  and  it  is  reformed  over 
Cu0/Zn0/Al203  catalyst  via  the  following  overall  reactions  [30]: 

CH3OH  +  H20^>C02  +  3H2.  (reforming)  (5) 


Some  portion  of  the  methanol  also  decomposes  directly  to  CO 
and  H2  via  the  reaction: 

CH3OH-^CO  +  2H2.  (decomposition)  (6) 

In  addition,  the  water-gas  shift  reaction  adjusts  the  composition 
of  the  product  gas: 

C0  +  H204^C02  +  H2.  (shift  reaction)  (7) 

It  is  documented  that  the  reforming  reactions  can  be  considered 
irreversible,  and  the  shift  reaction  could  be  neglected  without  a 
substantial  loss  in  accuracy  [30].  Therefore,  the  forward  reaction  of 
the  reforming  reaction  and  decomposition  reaction  are  under  con¬ 
sideration  in  this  paper.  The  chemical  reaction  rates  of  the  mixture 
on  the  catalyst  are  described  as 

rR  =  wcat/<RCi  (8a) 

rD  =  wcat/<D  (8b) 

where  wcat  is  the  catalyst  density  on  the  channel  wall,  C\  represents 
the  molar  concentration  of  methanol  near  the  catalyst  surface,  /<r 
and  kD  are  the  rate  constants  for  reforming  reaction  and  decompo¬ 
sition  reaction,  respectively.  The  reaction  rate  constants  depend  on 
the  properties  of  the  catalyst  and  reaction  condition.  The  quantities 
for  the  catalyst  F3-01  (Cu0/Zn0/Al203 )  can  be  related  to  Amphlett’s 
rate  constants  [30],  which  are  derived  form  an  Arrhenius  reaction: 

kR  =  CR[AR  +  BRln(SMR)]exp(-J|)  (9a) 

kD  =  CdAd  exp  (-j^r)  (9b) 

where  SMR  is  the  molar  ratio  of  steam  to  methanol,  Ar,  Br  and 
Ad  are  Amphlett’s  constants  for  reforming  and  decomposing  reac¬ 
tion,  respectively.  CR  and  CD  are  the  empirically  correction  factors 
for  reforming  and  decomposing  reactions  with  the  catalyst  F3- 
01(CuO/ZnO/Al2O3),  which  are  5.5  and  3.5,  respectively  [31].  The 
molar  rate  of  creation/destruction  of  species  i  can  be  expressed  as 


II 

1 

I 

o* 

(10a) 

& 

I 

II 

£ 

(10b) 

r3  =  3rR  +  2  rD 

(10c) 

** 

II 

3 

(lOd) 

r5  =  rR 

(lOe) 

where  the  subscripts  indicate:  l-CH3OH,  2-H20,  3-H2,  4-CO,  and 
5-C02. 
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Table  1 

Geometry  parameters  for  constructal  and  serpentine  microreactors. 
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Constructal 

k 

0 

1  2 

3 

4 

5 

6 

L/c /mm 

10.00 

7.07  5.00 

3.54 

2.50 

1.77 

1.25 

Serpentine 

Overall  dimensions 

wkl  mm 
h/ mm 

Total  flow  length/mm 

w/ mm 
h/ mm 

Total  channel  length/mm 

2 

1.12  0.72 

Lx  =  20  mm,  Ly  =  27  mm,  Lz 

0.50 

1 

62.3 

2 

1 

110 

=  3  mm 

0.36 

0.27 

0.20 

Fig.  3.  Computation  domain  of  the  microreator  with  constructal  tree-shaped  net¬ 
work. 


Considering  that  the  microreactor  with  constructal  tree-shaped 
network  is  symmetric  to  the  plane  ofy  =  0,  half  of  the  microreactor, 
shown  in  Fig.  3,  is  selected  as  the  computational  domain  to  sim¬ 
plify  the  calculation.  And  the  computational  domain  is  the  whole 
microreactor  for  the  serpentine  flow  pattern  reactor.  In  addition, 
the  following  assumptions  are  applied  to  the  model: 

1 .  Steady  state. 

2.  Negligible  gravity. 

3.  Mixture  of  the  steam  and  methanol  behaves  like  an  ideal  gas. 

As  to  the  flow  in  microchannel,  the  validity  of  the  continuum 
assumption  can  be  checked  by  evaluating  the  Knudsen  number 
(Kn): 

Kn=j  (11) 

where  X  is  the  molecular  mean  free  path  and  l  is  the  characteristic 
geometric  dimension  of  the  channel.  For  the  smallest  channel  size, 
i.e.  the  6th  hydraulic  diameter  d6  =  0.33  mm,  the  Knudsen  number 
are  still  far  lower  than  0.001 ,  which  assures  the  validity  of  the  con¬ 
tinuum  and  Navier-Stokes  equations  for  the  systems  considered  in 
this  work  [32].  In  addition,  the  maximum  Reynolds  numbers  under 
the  calculation  conditions  is  no  more  than  90  in  the  simulation,  so 
the  laminar  flow  of  the  mixture  is  considered. 

Thus,  the  governing  equations  for  the  mixture  in  the  microreac¬ 
tors  are: 

Mass  conservation 


V-(pmu)  =  0 


(12) 


Fig.  2.  Schematic  of  reactors:  (a)  reactor  with  constructal  tree-shaped  network;  (b) 
constructal  tree-shaped  network;  (c)  reactor  with  serpentine  flow  pattern. 


where  the  density  of  the  gas  mixture,  pm,  are  obtained  by  ideal  gas 
law. 

Momentum  conservation 


V  •  (Pmuu)  =  -Vp  +  V  ■  p,mVii 


(13) 
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where  the  viscosity  of  the  gas  mixture,  /xm,  is  computed  on  the  basis 
of  mass  fraction  average  of  the  pure  species 

Mm  =  Y;YilM  (14) 

i 

Species  conservation 

V-(pmuYi)  =  -V.Ji  +  Miri  (15) 

where  Yt  is  the  mass  fraction  of  species  i,  and  M,-  is  the  molecu¬ 
lar  weight  of  species  i.Ji  is  the  diffusion  flux  of  species  i,  which  is 
computed  by  solving  the  Maxwell-Stefan  equations  [33]. 

Energy  conversation 

V  .  (pmcpuT)  =  V  ■  kmVT  -  V  •  -  A HRrR  -  A HDrD  (16) 

j 

AHr,  A Hd  are  the  enthalpy  of  reforming  reaction  and  decomposi¬ 
tion  reaction  respectively  and  hj  is  the  enthalpy  of  species  j  (J  =  1,2, 
. . .,  5).  The  specific  heat  capacity  cp  and  thermal  conductivity  km  of 
the  gas  mixture  are  obtained  from  the  mass  fraction  average  of  the 
pure  species 


5 


y>p, 

i 

(17) 

5 

(18) 

1 


For  the  solid  part  of  the  microreactors,  there  is  no  mass  transport 
and  reaction  occurred.  The  energy  conversation 

V2T  =  0  (19) 


"in  (m/s) 


"in  (rnL/h) 


Fig.  4.  Comparisons  of  methanol  conversion  between  numerical  results  and  exper¬ 
imental  data:  (a)  tubular  channel  and  (b)  serpentine  flow  channel. 


Boundary  conditions  are  required  to  close  the  mathematical  for¬ 
mulation.  At  the  inlet  of  the  microreactor,  the  velocity,  temperature 
and  the  steam  to  methanol  molar  ratio  (SMR)  of  the  methanol  steam 
mixture  are  specified.  The  pressure  at  the  outlet  is  assumed  to  be 
atmospheric  pressure.  Adiabatic  boundary  conditions  are  applied 
to  the  reactor’s  external  surface  except  the  bottom  surface  which  is 
assumed  as  constant  wall  temperature.  Symmetry  boundary  condi¬ 
tions  are  applied  to  the  plane  aty  =  0,  which  is  the  symmetric  plane 
of  the  microreactor.  The  channel  wall  is  assumed  to  be  no-slip  and 
impermeable,  and  the  continuities  boundary  of  temperature  and 
heat  flux  is  applied  as  the  conjugate  boundary  condition  to  couple 
the  energy  equations  for  the  fluid  and  in  the  wall.  It  is  assumed  that 
there  is  no  catalyst  deactivation  in  the  simulation. 

3.2.  Numerical  solution 

The  set  of  governing  differential  equations  with  the  boundary 
conditions  described  above  for  the  microreactors  is  numerically 
solved  with  a  finite  volume  scheme,  and  the  pressure-velocity 
coupling  is  accomplished  by  using  the  SIMPLE  algorithm.  For  the 
complex  constructal  structure  evaluated,  a  structured  mesh  based 
on  hexahedron  grid  elements  is  applied  to  arrive  at  a  solution.  The 
grids  near  the  solid/fluid  interface  and  at  bifurcation  are  specially 
refined  with  consideration  of  the  effects  of  both  the  boundary  layer 
flow  and  bifurcation.  The  resulting  system  of  algebraic  equations 
is  solved  using  the  Gauss-Seidal  iterative  technique,  with  succes¬ 
sive  over-relaxation  to  improve  the  convergence  time.  A  careful 
grid-independent  study  is  performed  in  order  to  produce  grid- 
independent  results.  The  convergence  criteria  for  the  normalized 
residuals  for  each  variable  are  restricted  to  be  less  than  1 0-6  besides 
that  the  normalized  residuals  for  energy  is  restricted  to  be  less 
than  10-12.  In  addition,  the  solution  is  regarded  as  convergent,  not 
only  by  examining  residual  levels,  but  also  by  monitoring  relevant 
integrated  quantities  and  checking  for  heat  and  mass  balances. 


3.3.  Case  validation 

The  experiment  of  methanol  conversion  for  a  wall-coated 
tubular  microreactor  with  catalyst  thickness  of  100p,m  at  503  K 
conducted  by  Lee  et  al.  [13]  is  numerically  estimated  using  the 
above  model.  In  addition  to  the  tubular  microreactor,  the  present 
model  is  also  verified  by  a  miniaturized  microreactor  with  a  more 
complex  geometry,  serpentine  channel  [  1 0].  As  shown  in  Fig.  4,  the 
numerical  results  for  these  two  types  of  structure  are  in  agreement 
with  the  experimental  data  available  in  the  literature  [10,13].  This 
agreement  between  the  numerical  results  and  the  experimental 
data  verifies  the  present  model  is  reasonable. 

4.  Numerical  results  and  discussion 

In  order  to  evaluate  the  performance  of  the  reactor  for  methanol 
steam  reforming  reaction,  the  conversion  of  methanol  and  the  car¬ 
bon  monoxide  (CO)  concentration  in  the  product  (mixture  of  H2, 
C02  and  CO  at  the  outlet)  are  the  two  main  factors.  The  conversion 
rate  of  methanol  for  the  microreactors,  x,  is  calculated  by  following 
equation: 


where  C\  in  and  Ci)0Ut  are  the  inlet  and  outlet  molar  concentrations 
of  methanol,  and  the  subscript  1  refers  to  methanol. 

Fig.  5  presents  the  methanol  conversion  rate  and  CO  concentra¬ 
tion  as  a  function  of  steam  to  methanol  molar  ratio  (SMR)  in  the 
constructal  reactor  under  the  conditions  of  the  bottom  wall  tem¬ 
peratures  Tw  =  503  K,  gas  hourly  space  velocity  GHSV  =10,000  h-1 , 
and  that  the  feeding  mixtures  of  steam  and  methanol  is  of  the  same 
temperature  as  the  bottom  wall.  As  shown  in  the  figures,  a  larger 
SMR  not  only  improves  methanol  conversion  rate  but  also  reduces 
CO  concentration  in  the  product.  It  should  be  noted  that  the  benefit 
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Fig.  5.  Effect  of  molar  steam  methanol  ratios  on  the  performance  of  constructal 
reactor  (TW  =  503K  and  GHSV=1 0,000  h_1):  (a)  methanol  conversion  rate  and  (b) 
CO  concentration. 


is  offset  by  the  increased  heating  requirement  for  vaporizing  the 
additional  water  in  the  feed.  So  a  compromise  value  of  SMR=  1.1  is 
utilized  in  the  following  simulation. 

Fig.  6  plots  the  variation  of  the  species  mole  fractions  along  the 
flow  direction  in  tree-shaped  microreactor.  As  shown  in  the  fig¬ 
ures,  the  variation  range  of  species  mole  fractions  in  the  inlet  layer 
is  larger  than  that  in  the  outlet  layer.  In  other  words,  the  reaction 
rates  of  the  reforming  and  the  decomposition  are  faster  in  the  inlet 
layer  of  tree-shaped  network  than  that  in  the  outlet  layer  of  tree¬ 
shaped  network,  which  is  induced  by  larger  mole  concentration  of 
methanol  in  upstream.  In  addition,  in  the  inlet  layer  of  tree-shaped 
network,  the  variation  of  species  mole  fractions  mainly  occurs  in 
the  region  from  the  3th  branch  to  the  6th  branch.  This  fact  implies 
that,  in  the  tree-shaped  architecture,  with  the  increasing  branches, 
the  factors,  including  the  reduction  of  reactant  velocity  and  larger 
surface-to-volume  ratio  owing  to  the  reduction  of  channel  size, 
lead  to  enhance  the  methanol  steam  reforming  process.  In  order  to 
present  a  clear  understanding  of  reforming  process  in  tree-shape 
network,  Fig.  7  illustrates  the  mole  fraction  of  hydrogen  at  typical 
cross  sections  including  both  in  the  inlet  and  outlet  layers  of  tree- 
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Fig.  6.  The  variation  of  mole  fractions  along  the  flow  direction  in  constructal 
microreactor  ( Tw  =  503  K,  GHSV  =  1 0,000 h-1,  SMR=  1.1). 


Fig.  7.  Distributions  of  H2  mole  fraction  (the  same  condition  as  Fig.  6):  (a)  atz  =  0  mm 
(inlet  layer)  (b)  at  z  =  1 .3  mm  (outlet  layer). 


shaped  network.  It  is  indicated  by  Fig.  7  that  the  variation  of  mole 
fractions  of  H2  mainly  occurs  in  the  region  of  latter  branch  level. 

A  comparison  between  the  conversion  in  microreactor  with 
tree-shaped  network  and  microreactor  with  serpentine  flow  pat¬ 
tern  with  respect  to  the  gas  hourly  space  velocities  and  bottom  wall 
temperatures  is  presented  in  Fig.  8.  As  seen  from  the  figure,  the 
methanol  conversion  in  the  constructal  microreactor  is  more  than 
10%  over  that  of  serpentine  reactor  at  the  same  gas  hourly  space 
velocity  and  bottom  wall  temperature.  The  superiority  is  driven  by 
the  fact  that,  the  tree-shaped  configuration  provides  reaction  space 
with  larger  surface-to-volume  ratio  owing  to  the  reduction  of  chan¬ 
nel  size  in  the  latter  branches,  and  the  reactant  velocity  reduces 
in  the  latter  branches.  It  is  also  indicated  by  the  figure  that,  both 
conversion  rate  in  the  two  types  of  reactors  decrease  as  the  space 
velocity  increases,  since  the  methanol  passes  through  the  reactors 
quickly  and  has  less  time  to  contact  with  the  catalyst  with  a  larger 
space  velocity.  And  the  conversion  in  these  two  type  reactor  obvi¬ 
ously  increases  with  the  increase  of  the  bottom  wall  temperature. 
Therefore,  a  higher  methanol  conversion  can  be  obtained  under  a 
small  space  velocity  or  a  higher  reaction  temperature. 

The  CO  concentrations  in  both  the  constructal  tree-shaped  net¬ 
work  microreactor  and  serpentine  microreactor  with  respect  to  gas 
hourly  space  velocities  and  the  bottom  wall  temperatures  are  plot¬ 
ted  in  Fig.  9.  As  shown  in  the  figure,  both  the  CO  concentrations 
in  the  product  of  the  two  types  of  microreators  decrease  with  the 
increasing  space  velocity  but  increase  with  the  increasing  reaction 
temperature,  which  are  of  the  same  trend  as  the  conversion  rate. 
In  other  words,  a  higher  methanol  conversion  rate  is  accompanied 
with  a  higher  CO  concentration  in  the  product.  And  the  CO  con¬ 
centration  in  the  product  of  the  constructal  microreactor  is  a  little 
higher  than  that  of  serpentine  microreactor. 

Fig.  10  compares  the  total  pressure  drop  between  the  construc¬ 
tal  tree-shaped  and  serpentine  reactors.  In  the  microreactor  with 
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Fig.  8.  Comparison  of  the  conversions  rate  between  constructal  and  serpentine 
microreactors:  (a)  Tw  =  503  K,  (b)  GHSV  =  10,000  hr1 . 


Fig.  10.  Pressure  drop  versus  various  gas  hourly  velocities  at  Tw  =  503  K. 

tree-shaped  network,  the  local  pressure  loss  is  induced  by  a  large 
number  of  the  bifurcations,  and  the  reduction  of  channel  size  in 
the  latter  branches  also  cause  pressure  losses.  However,  as  shown 
in  Fig.  10,  the  pressure  drop  in  the  microreactor  with  construc¬ 
tal  tree-shaped  network  is  not  greater  than  the  microreactor  with 
serpentine  flow  pattern  at  various  gas  hourly  space  velocities.  The 
primary  reason  is  that  the  constructal  tree-shaped  network  acts  as 
an  optimized  fluid  transport  path  which  can  effectively  reduce  the 
flow  distance.  For  the  constructal  reactor  considered  in  this  paper, 
the  flow  distance  is  just  56.6%  of  the  corresponding  serpentine 
reactor.  The  reduction  of  flow  distance  makes  the  tree-shaped  reac¬ 
tor  still  possess  almost  the  same  pressure  drop  as  the  serpentine 
reactor. 

5.  Conclusion 


GHSV  (h'1) 


Fig.  9.  Comparison  of  CO  concentrations  in  the  product  between  the  constructal 
and  serpertine  reactors:  (a)  Tw  =  503  K  and  (b)  GHSV=  10,000  h_1. 


The  constructal  tree-shaped  network  is  introduced  into  the 
design  of  the  methanol  steam  microreactors  in  the  context  of  opti¬ 
mization  of  flow  configuration.  A  three-dimensional  model  for 
methanol  steam  reaction  in  this  new  designed  microreactor  is 
developed  and  numerically  analyzed.  The  methanol  conversion, 
CO  concentration  in  the  product  and  the  total  pressure  drop  of 
the  gases  in  the  reactor  with  constructal  tree-shaped  network  are 
evaluated  and  compared  with  those  in  a  reactor  with  serpentine 
flow  pattern,  which  is  of  the  same  reaction  volume  and  overall 
dimensions.  The  conclusions  can  be  summarized  as  follows: 

(1)  Compared  with  the  serpentine  reactor,  the  constructal  reac¬ 
tor  possesses  a  higher  methanol  conversion  rate  accompanied 
with  a  higher  CO  concentration  in  the  product.  The  methanol 
conversion  of  the  constructal  microreactor  is  more  than  10% 
over  that  of  serpentine  reactor.  This  superiority  is  driven  by 
the  fact  that,  the  tree-shaped  reactor  configuration,  which  acts 
an  optimizer  for  the  reactant  distribution,  provides  a  reaction 
space  with  larger  surface-to-volume  ratio  and  the  reduction  of 
reactant  velocities  owing  to  the  reduction  of  channel  size  in  the 
branches. 

(2)  In  the  tree-shaped  network,  despite  the  fact  that  the  bifur¬ 
cations  induce  extra  local  pressure  loss,  and  the  reduction  of 
channel  size  in  the  latter  branching  level  also  causes  pressure 
losses,  however,  the  reduction  of  flow  distance  makes  the  reac¬ 
tor  with  constructal  tree-shaped  network  still  possess  almost 
the  same  pressure  drop  as  the  serpentine  reactor. 
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